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Detection of a weak magnetic field via cavity enhanced Faraday rotation 
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We study the sensitive detection of a weak static magnetic field via Faraday rotation induced 
by an ensemble of spins in a bimodal degenerate microwave cavity. We determine the limit of the 
resolution for the sensitivity of the magnetometry achieved using either single-photon or multiphoton 
inputs. For the case of a microwave cavity containing an ensemble of Nitrogen-vacancy defects in 
diamond, we obtain a magnetometry sensitivity exceeding 0.5 nT/\/Hz, utilizing a single photon 
probe field, while for a multiphoton input we achieve a sensitivity about lfT/\/Hz, using a coherent 
probe microwave field with power of Pin = 1 nW. 

PACS numbers: 33.57.+C, 07.55.Ge, 42.50.Pq, 85.70.Sq 


I. INTRODUCTION 

Sensing of magnetic fields with an extremely high sen¬ 
sitivity has attractive applications in various areas of 
science and technology m- Magnetic fields can be 
measured using a variety of techniques including a su¬ 
perconducting quantum interference device (SQUID) [T], 
magneto-materials [1], atoms [ME] and color defect cen¬ 
tres in diamond mn]. The SQUID magnetometers 
achieve record sensitivities [1] , but require extremely low 
temperatures to maintain superconductivity [T]. Re¬ 
cently, atomic magnetometers have demonstrated sub- 
femtotesla sensitivity, approaching the record sensitivity 
of SQUID sensors [BHH]. Diamond magnetometers ex¬ 
ploiting nitrogen vacancy (NV) defect centres in diamond 
offers detection of magnetic field signals both with high 
spatial accuracy [Mis] as well as high field sensitivity 
down to sub-pT/vTjz [3 [TT] . In this work, we describe a 
magnetometer approaching sub-fT/V^Hz sensitivity with 
nW-order input power that exploits the Faraday rotation 
of microwave (mw) photons induced by a static magnetic 
held. 

In the Faraday effect the polarization of electromag¬ 
netic helds traveling in magneto-optical material can be 
rotated by applying a static magnetic held. This Faraday 
rotation can also be well explained by a quantum mecha¬ 
nism |16j . Based on this quantum understanding Faraday 
rotation has been proposed to detect quantum huctua- 
tions m and induce giant phase modulation |I8j . Giant 
optical Faraday rotation has also been observed mm- 
The polarization rotation due to the Faraday effect has 
been well studied for classical optical magnetometers but 
a corresponding discussion of the sensitivity limit is so far 
absent in the literature. 

In this paper we investigate the sensitivity limit of a 
microwave version of a Faraday magnetometer. In our 
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scheme the polarization of microwave probe photons is 
rotated by an ensemble of spins coupled to a microwave 
cavity. By measuring the subsequent rotation of the pho¬ 
ton’s polarization we are able to detect a weak static 
magnetic held applied to the spins. Thanks to the low 
frequency of microwave photons and the lower vacuum 
energy huctuations, as compared with the optical pho¬ 
tons, we hnd that the microwave rehectance Faraday 
magnetometer is ultrasensitive in comparison with op¬ 
tical magnetometers, given the same probe held power. 


II. SYSTEM AND MODEL 

Before describing the details of our system and model 
we hrst present the main idea of our scheme to measure 
a weak magnetic held. Our scheme involves an ensemble 
of spins coupled to a mw cavity. A horizontal-polarized 
probe mw held is incident to the cavity. This mw input 
held can be decomposited into the cr+- (right circular- 
polarized) and cr_-polarized (left circular-polarized) com¬ 
ponents. The mw cavity also possesses ct_|_- and cr_- 
polarized modes which are degenerate in frequency. The 
two orthogonally polarized cavity modes couple to two 
separate transitions of spins contained in the cavity. 
When a weak static magnetic held is applied to the en¬ 
semble of spins it shifts the two transition energies up 
or down oppositely and subsequently causes diherent de¬ 
tunings between the transitions and the cavity modes. 
Thus, the rehected helds suffers different phase shifts 
yielding a Faraday rotation angle. As a result, in the 
output the superposition held of the two polarization- 
orthogonal modes includes both horizontal and vertical 
polarized photons. Measuring these photons we can es¬ 
timate the magnetic held with a high precision. 

We now describe in more detail our scheme for mag¬ 
netic held sensing. In the scheme an ensemble of NV cen¬ 
ters couples to the mw cavity as shown in Fig. We as¬ 
sume that the mw cavity supports right and left circular- 
polarized ((7+ and CT_-polarized) cavity modes, {a+,a_}, 
and these modes are degenerate with resonance frequency 
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FIG. 1. (Color online) Schematics of setup for detection of a 
weak static magnetic field B. An ensemble of spins couples 
to a microwave (mw) cavity. The horizontal polarized mw 
field ain inputs into the mw cavity and then is reflected off 
the cavity (Oout) to the detector. The output field reflected 
off the cavity is detected by the polarization resolving photon 
detector D. The polarization of the output field is rotated by 
an angle ip_F due to the Faraday effect. The cavity has the 
resonance frequency uJr and the intrinsic loss Ki). The input 
mw field couples to the cavity with an external coupling rate 

^ex • 

uir, and intrinsic loss kj. A horizontal (H)-polarized co¬ 
herent mw field Qfin with frequency Win enters the cavity 
with an external coupling rate of Kex and then yields an 
output field aout- This input field, a-m, can be decom¬ 
posed into cr+ and cr_-polarized components driving the 
cavity modes. The (T+ and (T_-polarized cavity fields suf¬ 
fer different phase shift because they couple to different 
transition within the spin ensemble (Fig. [^. Due to the 
interaction between the spins and the mw fields in cav¬ 
ity, the polarization of the output field Oout is rotated by 
the Faraday angle, (pp, and hence aout includes a vertical 
(V)-polarized component. For a clear discussion of the 
relation between the input and output amplitudes and 
polarizations we define the annihilation and creation op¬ 
erators, c and ifi, for a quantized radiation field. We use 
the notation mull], to denote a y—polarized field 
of mode <f, where x is the unit polarization vector. The 
polarization vector x can be H, V, (t+, crT or an arbitrary 
linear polarization 6 = cos0H -|- sin0V with 0 G (—7r,7r]. 
Next we will discuss the magnetometer sensitivity in the 
case of single-photon and multiphoton measurements. 

In the case of a very weak input a single-photon probe 
mw field incidents to the cavity and we consider all pos¬ 
sible outputs. The output mw field is no longer perfectly 
horizontal polarized due to the Faraday effect and a po¬ 



FIG. 2. (Color online) Configuration describing the interac¬ 
tion between the cavity mode a± and an ensemble of NV“ 
centers. Two degenerate cavity modes, the right circular- 
polarized mode a+ and the left circular-polarized mode a_, 
are detuned from the zero-strain splitting by Aq. We shift the 
transition \ms = 1) -f-)- \ms = 0)(|ms = —1) \ms = 0)) by a 
bias static magnetic field Bq, which creates a static frequency 
shift A = fiBgeBo. The weak static magnetic field SB causes 
another small shift S — ^BQeSB. 

larization sensitive mw photon detector connected to the 
output port may detect a horizontal or vertical polarized 
photon or a “dark count” with associated probabilities. 
The “dark count” implies that the input photon is lost 
to the environment or is absorbed by the materials. The 
limit of the sensitivity can be determined by the Fisher 
information. In classical measurement, we input a weak 
coherent probe mw field with power Pin into the cavity 
and only measure the intensity of the vertical-polarized 
output component /^. The horizontal-polarized compo¬ 
nent is filtered from the output. 

Our scheme involves the interaction of mw modes in 
the cavity and an ensemble of spins, taken here, for 
example to be NV centers. In Fig. we graph the 
ground state triplet of the NV defect in the presence of 
a static bias field. The mw cavity modes a+ and a_ 
drive the magnetic transitions \ms = 1) O |tos = 0) 
and \ms = —1) O \ms = 0), respectively. In the ab¬ 
sence of any magnetic field, the zero-strain splitting of 
the NV centers is assumed to be P « 2.78 GHz. Un¬ 
der the bias magnetic field Bq and the signal field to be 
sensed, SB, the levels of \ms = 1) and |tos = —1) are 
shifted up or down by A -|- <5, respectively. The magnetic 
field is applied along the 0 -axis of the spin crystal, and 
results in (5 = and A = ^bq^Bq, where /i^ = 14 

MHz • mT“^ is the Bohr magneton and « 2 the elec¬ 
tron spin g-factor. The interaction Hamiltonian describ¬ 
ing the coupling between the spins and cavity takes the 
form 

Hi='^ (gi,i?a+|TOs = 0 )j(tos = 1| + H.c)j 

3 

+ ^ {gj^La3_\ms = Q)j{ms = -1| -|- , 

3 

where {9j,L) is the coupling rate between the cavity 
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mode a+ (a_) and the transition \ms = 1) O \ms = 0) 
(I ms = —1) -O' I ms = 0)) of the jth spin. For sim¬ 
plicity we assume that = 9 for all spins. 

For an ensemble of N spins we can apply the Holstein- 
Primakoff transformation [53] to define the collective 
operators c+ = \'^s = = 1| and 

c_ = = —1|, which allows one 

to consider a large number of spins as a generalized har¬ 
monic oscillator coupled to the cavity mode with a collec¬ 
tively enhanced rate G = VNg- The Holstein-Primakoff 
transformation requires that the spin ensemble is highly 
polarized in the |ms) state which can be achieved via op¬ 
tical pumping. Hence, the interaction Hamiltonian now 
becomes 

Hi = (Galc+ + G*alc- + H.c^ . (1) 

The free evolution Hamiltonian of spins under the mag¬ 
netic fields Bq and 5B is given by 

^^spin = j -b S)Szj , 

j j 

where Szj is the z component of the spin-1 operator, and 
Dj the zero-strain splitting of the jth spin. If we neglect 
the inhomogeneous broadening of Dj due to different lo¬ 
cal strains and assume Dj = D, then we have 

-ffspin = D{c\c+ + cLc-) -b -b ^)(c^c+ - cLc-). (2) 


the input and output in the H and V-polarized basis is 
governed by a scattering matrix as 




= Sr 


dPU 


(5) 


with 
Sr = 


HH 

^VH 


ir 


HV 

vv 


r++r- 

2 

. —r_ 




’~+ 




( 6 ) 




Specifically, when a H-polarized mw probe field = 
= ag/V2(cT+ -b CT_) is input into the cav¬ 
ity, the output field takes the form doutXout = + 

r_)H — i^(r+ — r-)V. The corresponding output 
power spectrum of the vertical and horizontal polar¬ 
ization are given by S'^(win) = i?v(a“^(-Win)dS(‘^in)) 
and Sjj(u;in) = i?H(dS^(-Win)dS(‘^m)) [2Q], respectively, 
where By = \ryfj\^{RYi = We define the re¬ 
flection r± = f = gj. = Ifbtizhbd ^ 

The polarization of the output field is rotated by Faraday 
angle Lpp = with respect to the H-polarized input 

field. For a V-polarized input field a^V = m^/-\/2((T+ — 
(T_), the output field is f%(r+ — r-)H -b ^(?'+ -br_)V. 


III. LIMIT OF MAGNETOMETRY 
SENSITIVITY 


We consider input probe fields d^, each with frequency 
Win. The full Hamiltonians for (T+ and CT--polarized cav¬ 
ity modes and spin ensemble takes the uniform form 

=Wnd±a± -b D(^j.c± ± (A -b (5)c^c± (3) 

+ {Gaic± + G*cia±) 

+ iV2Kex(d^d± - , 

where the tT+((T_)-polarized input uses the -b(—) sign, 
and 7 is the decoherence rate of the spin ensemble for 
each transition. Using these Hamiltonians and the input- 
output relation |53J , we can calculate the amplitude 

of the reflection in the frequency domain at w = Win to 
be 

r±(wi„) = -1 -b I I G" ’ 

where the detunings Ar = ujr — Win and Aq = D — Win. 
Win is the carrier frequency of the input probe field of 
Qfin. Equation Q is valid for both a single-photon and 
classical probe field when the number of photons in the 
cavity is much smaller than the number of spins [2BH28j . 
We convert the input and output fields from the <7+ and 
CT_-polarized basis to the H and V-polarized basis by the 
relations a+ = {H — — {H + iV)ly/2 [or 

H = ((7+ -b n_)/-\/2, V = f(i7+ — a-)l\/2] [531 . Thus, 


A. Single-photon input 

If we assume that our detector has a unit quantum effi¬ 
ciency ry = I and input a single photon into the setup, the 
probability to detect an output photon polarized along V 
or H can be determined as = 5'^/S'in or = Sfj/S-m 
with 5'in(w) = (a;jj(—Win)ain(win)) at w = Win- It is also 
possible that no photon clicks the detector due to the 
“dark count” indicating the loss of the photon before 
clicking the detector. This latter probability due to the 
loss of photon from cavities is P 0 = 1 — P^ — Pg. The 
probabilities for three possible outputs are 

P- =f‘^sin^{(pF) + Sr^cos^{<fF) = , 

Pg = f^ cos^{(f f) + Sr"^ = Irggl"^ , (7) 

P 0 = 1 - (f^ -b (5r^) . 

We are interested in the sensitivity of the system to 
perform a measurement of the static signal field 6B, 
, which means how fast the probabilities to de¬ 
tect a photon in state ^ = {V, H, 0} change for a certain 
signal SB. We first examine the limit of the magnetom- 
etry sensitivity using a single-photon input mw pulse. 
Consider a single-photon probe prepared in an initial 
quantum state /9(0) that this evolves to a state p{Tm) 
when exposed to the signal 5 = after a measure¬ 

ment time Tm- This generates three possible outcomes 
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with probabilities with ^ G {i?, F, 0}. To evaluate 
the performance of sensitivity of our setup, we rescale all 
parameters by Hi. Generally, the maximum amount of 
information about 5 that can be extracted from the po¬ 
larization dependent measurement is given by the Fisher 
information [32] 


Fi{5B) = 


f ^^B9t 


V 


E 


1 fdP{i\5/Ki) 


i=V,H,(Z> 


P{£.\6/Hi) V d5/K, 


= {9-B9ef ^ 


1 






pm \ ds j 


This leads to the Cramer-Rao bound 
AB > ^ 


\/vmFi{5B) ’ 


( 8 ) 


(9) 


where Vm is the number of times the measurement is 
repeated. If the total measurement time is Ttotai and 
each measurement takes a time Tm, then we have the 
following sensitivity for a single-photon input 


ABspV^ > . (10) 

Fi{dBj 

In practice the time is determined by the full width at 
half maximum (FWHM) of the Fisher information when 
evaluated as a function of the signal SB (See Fig. 
below). 


B. Multiphoton input 


Generally, taken estimating the value of a parameter ip 
contained in a measurement operator M, the root mean 
square (RMS), A(^, of variance as an imprecision in the 
estimation can be obtained from [MH55] 


A(/3^ 


AM2 

]d{M)/d^\^ ’ 


( 11 ) 


where the variance AM^ = (M^) — (M)^, with the expec¬ 
tation values taken as the appropriate input state. We 
only detect the V-polarized output photons so we have 

M = aZm_ 

We now estimate the limit of the measurement sen¬ 
sitivity for the multiphoton input where the probe field 
is a weak coherent mw pulse with an input power Pjn 
corresponding to a mean photon number hin, the quan¬ 
tum expectation of the operator fiin = We have 

(dii/db) = 2^6x11111 = [371 [3H1 (din = llSJ)j 


where I/tex = 2/Cex is the photon decay rate into the 
associated outgoing modes m- Here we replace Tex 
with the duration of the probe pulse when calcu¬ 
lating riin m because the duration of probe pulse is 


limited by the bandwidth determined by Fisher infor¬ 
mation now. The spectrum of the input probe mw field 
*S'in(win) = (dJ^l(-Wi„)d(),(a;in)) = 2Kexni„ [501 IMj. We 
now focus on the vertical polarized output. 

In our setup, the input and output ports, and the cav¬ 
ity support both H- and V-polarized fields. As a result, 
both the H- and V-polarized noise can enter the input- 
output port and the cavity through the external cou¬ 
pling channel or the intrinsic lossy channel, and then are 
reflected to the detector. The operators denoting the 
noise entering the cavities are = 'CE(d'-i- + 
and IgV = i|g((T+ — due to noise entering 

the cavity via the external environment through the 
Kex channel, and = Ci^(ct+ -I- a-)l^/2 and = 
~^-)l due to internal loss channels within the 
cavities. We define nominal transmissions t± = 1 + r± 
for use below. 

The quantum Langevin equations for the cavity modes 
and spin operators c± now take the form 

d± ={—iu)r — K)a± — iGc± + (I2a) 

A ^ex^E P- 7 

c± =[—i{D ± A ± (5) — 7/2]c± — iGa± , (12b) 

where k = Hex + Ki is the total decay rate of the cavity. 
The modes c± are Holstein-Primakoff transform of the 
collective spin operators, and we don’t consider noise en¬ 
tering them. The internal noise at w = Win, and 
also enter the output ports through the scattering matrix 



so that 


aoutF 

al.V 


= St 




The output fields are 


oMFI 
alV 




(14) 


^out — (“ki + Ie A *Ce ) 



The vertical polarized output includes four contributions: 
(I) the input probe field, fib; (II) the external H-polarized 
noise (III) the internal H-polarized noise (IV) 
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the external V-polarized noise (V) the internal V- 
polarized noise The total V-polarized output field is 
given by 


®out “ 



(16) 


Note that (r+ — r_)/2 = {t+ — t.)/2 and (r_|_ -|- r_)/2 = 
— 1 -l- (t+ -l- t.)l2. 

The whole device is put in an environment with a 
uniform temperature T. The cavity thermal noise op¬ 
erators are defined as 

with X G {H,V}. Their quantum expectations are 
(hg) = 2Kexhth and (fif) = 2Kihth [10], respectively, with 
fi-th = xji^^rlKBT _ jg Boltzman constant. 

Note that only near-resonant noise enters the cavity and 
detector. Their variances are = (2Kex)^ (h-th + ^th) 
and = (2Ki)^ (h^j^ -I- hth) HO]- We also have (hin) = 
2Kexhin and Ah?j = (2Kex)^hin for the coherent input 
field. The quantum expectations and the variances of 
the H- and V-polarized thermal fields are equal. 

The small magnetic field SB to be measured is con¬ 
tained in the measurement operator M = The 

V-polarized output field is a noisy coherent state with 
the total thermal noise ^tot = + 

{’’"hh^e + (1 + f°P ^ coherent state 

I - We define = (ILItot) and have gD] 


= 2KexT’yrith 


1-H 



2 ' 


T 2/^exT*^vrth 


l-f-^ 


2^[r 


HHi 


H 



where 3?[a;] gives the real part of a complex number x. 
The quantum expectation of the operator M reads 


(M) = 2Kex^’y”m + , (18) 

where Thus, the variance of the output 

can be evaluated as 


AVT^ — (2Kex^y)^*^in(2n^/2KexT\/ -f 1) -l- (n| -f n^) . 


When the mean photon number of thermal noise entering 
the detector is much larger than unity, and the photons 
in the probe field is much larger than this thermal noise, 
i.e. hin ^ hj ^ 1, to a good approximation, we have 


AM^ = AKey^PyTlin (uj + Ke^Py) , (19) 


and also 


[ dS J 


(d^ 

\ ds ' ■ 


^ ' nl 


( 20 ) 


Substituting Eqs. (19) and (20) into Eq. ( jlTj ), and 
using S = ^sgeSB, 2Kex7ex = 1, we obtain the sensitivity 
limit for the multiphoton input 


ASmpV^t^ > 


(2n^ + 2Kex^y 


> 


V^^^ex^in 

1 y^CthKsT + P^hlOr 






( QP 'N ^ 

is the nom¬ 
inal Fisher information of the vertical-polarized 

2" 


output field and Cth = 2P, 


V 


2P, 


H 


1 


l+P^+2U[rgg] 


i^y 


!+(*)■ 


Here we applied 


the relation rith^r ~ KbT for nth ^ 1. In the case of 
Ki <C Kex, the limit becomes 


Ai?MP-\/'rtot ^ 


1 




)KqT + p^hujj- 


\/ -Pl.v(^) 


VK 


( 21 ) 

The multiphoton limit of Eqs. (21) immediately shows 
us that the sensitivity is proportional to \/T/Pin for 
KbT 3> hlJJr^ and can be improved by cooling the res¬ 
onator or increasing the probe power. 

It is very interesting to note that the sensitivity limit 
in the multiphoton measurement Eq. (21), is inversely 
proportional to \/Pin oc \/Tml hin ■ It means that the 
sensitivity cannot be improved by increasing the dura¬ 
tion of the probe pulse when keeping the input power 
constant. 

The multiphoton limit also shows an important ad¬ 
vantage for sensing using a mw cavity with an input 
mw field of frequency Wmw over an optical system with 
Wo ^ Wmw Given the same Fisher information, the in¬ 
put power Pin and the environmental temperature we 
have IiWniw/2 <C KbT <C fvjJo/2, the sensitivity can be 

improved by ^ using the mw system in com¬ 

parison with the optical system. For typical parameters 
T = 70 K and Wq = 1.78 x 10^® rad/s corresponding to 
the wavelength A = 1064 nm, the improvement factor 
is about 10 and can be up to two orders as P —>■ 0 for 
Wmw = 27r X 3 GHz. 

Our system is too complex to provide an analytical 
form for the Fisher information and in what follows we 
calculate it numerically. 


IV. NUMERICAL RESULTS 

Figure shows the probabilities, Eqs. ([^, and their 
derivatives for the three possible outcomes (horizontal, 
vertical and no photon). The conditional probabilities 
P(^|i5) (^ = {V, P, 0}), are found to be symmetric 
with respect to 5 = 0. The probability P{0\6) has a 
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FIG. 3. (Color online) Probabilities P{^\5/Ki) (a) and deriva¬ 
tion of probabilities dP{^\S)/d5 (b) of three outputs as a func¬ 
tion of the magnetic-field induced level shift 5. Ab = = 

Aq = 0, A — 0,Kex — Ki,G — K.iy'f = 10“^. Blue lines in¬ 
dicates the probability of the vertical polarization (Py), red 
lines for the horizontal polarized outcome (Ph) and green lines 
for the probability of detecting zero photon (P 0 ). 


Lorentzian profile with a dip P{0\S) = 0 at (5 = 0. 
When (5 ~ 0, the probability to detect a horizontal po¬ 
larized photon is nearly unity. When the cavity and the 
spins are off resonant the probability to detect a verti¬ 
cal polarized photon increases and reaches the maximum 
P{V\6) = 0.25 at |i5| = OAdAm, and then decreases again. 
This nonzero value of 6 indicates that there is an optimal 
bias with A = S. Throughout the description below, we 
will replace the nonzero S with the bias A. 

The derivatives of these probabilities are antisymmet¬ 
ric with respect to <5 = 0. They reach abolsute maxima 
at similar positions, S ~ 0.2Ki. Therefore, when the spins 
are biased at A ^ 0.2Ati, our magnetometer is most sen¬ 
sitive to the weak magnetic signal. 

The Fisher information is crucially dependent on the 
conditional probability derivatives. At (5 = 0, the Fisher 
information has a deep drop, while it has a peak of 

Fi{S/Ki) = 29 at |A| = 0.07k,. The full FWHM 


defining a bandwidth for measurement is about O.dKi. 
Note that the resolution of a probing system is governed 
by the product of the average Fisher information and the 
bandwidth l/vm, see Eq. (10). This gives a sensitivity 
of AB-^Ttotai ^ 1-2 X 10“^\/Ki for G = Kex = Ki and 

7 = lO-^Kj. 

Now we find the optimal sensitivities. As shown in 
Fig.i we calculate the sensitivity for the single-photon 
input as a function of the external coupling Kex for dif¬ 
ferent spin-cavity coupling G. It can be seen that the 


FIG. 4. (Golor online) Fisher information as a function of 
the magnetic-field induced level shift S. Ar = Aq = 0, A = 
0, G = Kex = Ki,7 = 10~^Ki. 



FIG. 5. (Golor online) Sensitivity scaled by 
as a function of the external coupling Kex and the 
cavity-spin coupling G. G/ki = (0.02,0.1,0.2) = 

(solid blue line, dashed red line, dot-dashed yellow line). 

Ar = Aq = 0. 


optimal sensitivity of < 0.03 ^ ‘ can be obtained when 
G/k, ^ 0.1 and Kex/«^i > 7.For a large coupling, e.g. 
G /= 1, the sensitivity is low, > 0.1 ^^^ . A very small 
coupling, e.g. G/ki = 0.02, the sensitivity becomes lower. 
For mediate coupling rates, G/ki ~ 0.1, the sensitivity 
rapidly decreases to a limit of 0.027^^- at Hex /~ 10. 

To find the optimal spin-cavity coupling G, we calcu¬ 
late the Fisher information, shown in Fig. Clearly, 
the Fisher information Fj is large when 0.06 < G/k, < 
0.1, whereas the width increases from 6.4 x 10“^Ki to 
9.6 X 10“^Ki. This indicates that the sensitivity is opti¬ 
mal when G « O.Iki. 

For the multiphoton input, we are interested in the 
maximal Fisher information of the vertical-polarised out¬ 
put. Once known, we can estimate the sensitivity for any 
environmental temperature and any input power. Figure 
[ 7 ] shows the nominal Fisher information of the vertical- 
polarised output as a function of the detuning <5. The 

maximum value is about Fiy « 10® ^ ) yielding 

1.92 X lO-^V^ for Pin = 1 nW and T = 70 K. Note 
that the FWHM is small, about 4 x 10“®Ki. Thus, to 
achieve the highest sensitivity, the probe pulse need be 
long. 
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FIG. 6. (Color online) Fisher information scaled by ^ j 
as a function of the spin-cavity coupling G and the magnetic- 
field induced level shift S. Ar = A, = 0 and Kex = IOki. 
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FIG. 7. (Color online) Fisher information as a function of the 
intercavity coupling J and the detnning <5. A;, = Aa = Aq = 

0, 4 = 0, Kex/Ki = 10, 7/Ki = 10“^, G/Ki = 0.1. 



V. DISCUSSION OF EXPERIMENTAL 
IMPLEMENTATION 

The experimental implementation of our magnetome¬ 
ter crucially relies on the realization of a microwave 
cavity supporting both ct+- and (T_-polarized (or H- 
and V-polarized) modes, which are degenerate in fre¬ 
quency. Such microwave cavities require particular de¬ 
signs and have been demonstrated using transmission line 
resonators [4 11143 j and three-dimensional superconduct¬ 
ing Fabry-Perot microwave cavities [441l46j . The Q fac¬ 
tor of the transmission line resonator can be 100 — 10® 
[HUS], while that of the superconducting FP mw cavity 
can reach 10^^ [44]. Recently, Zhang et al. developed 
a photonic crystal(PC) cavity with a small mode vol¬ 
ume and a high Q (up to 10®), working at microwave 


frequencies [35] • This PC cavity can also possess two 
degenerate modes with orthogonal polarizations [49ll51j . 
In recent works gDia, the magnetic transition between 
the triplet ground states of NV centers has been selec¬ 
tively driven with circularly-polarized microwave fields 
in microstrip mw resonators. Thus many of the essential 
components for our scheme have already been experimen¬ 
tally demonstrated. 

Now we estimate the optimal sensitivity for detecting 
a weak magnetic held using an ensemble of NV centers 
coupled to a microwave cavity. Since the magnetic tran¬ 
sition frequency of NV is about 27r x 2.8 GHz, we require 
a microwave cavity with w,. = 27r x 2.8 GHz. The decay 
rate 7 of Im^ = ±1) of the ground state triplet has been 
measured to vary from a few MHz |52l [55] to O.OlHz [S31- 
[55] . Our setup with practical parameters 7 = 10“®Ki, 
Kex = lO^i and G = O.lKi, using a low-Q factor cross mi¬ 
crostrip resonator with Q = 100 [311-031; yields Kil2TT = 
28 MHz, and can provide a sensitivity of AHspV'Totai ^ 
5.2 nT/-\/Hz for a single-photon probe held. If an input 
power of Pin = 1 nW is applied and the environmental 
temperature is hxed to T = 70 K, the sensitivity can 
approach APmp V^totai ^ 1 fT/-\/Hz. Note that sen¬ 
sitivities scales as ^/Ki■ Therefore, if superconducting 
Fabry-Perot cavities with a mediate Q = 10® [33], corre¬ 
sponding to Ki/27r = 28 kHz, is applied, the sensitivities 
can be improved to APmpVA otai ^ 32.2 aT/vldz when 
Pin = 1 nW. 


VI. CONCLUSION 

We have studied the polarization rotation of linear- 
polarized microwave photons input into microwave cavi¬ 
ties coupled to an ensemble of spins. Measuring the Fara¬ 
day rotation of the output photons provides a method to 
ultrasensitively measure the strength of magnetic helds. 
The sensitivity limit of this microwave magnetometer is 
presented for both single-photon and multiphoton probe 
helds. The sensitivity of the magnetic held in the multi¬ 
photon measurement with an input power of Pin = 1 nW 
can be tens of aT/v^Hz. 
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